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We present a spectro-photometric catalogue of the evohition of supernova 1988Z 
which combines new and pubhshed observations in the radio, optical and X-ray bands, 
with the aim of offering a comprehensive view of the evolution of this object and 
deriving the total energy radiated since discovery. The major contribution to the 
total radiated energy comes at optical to X-ray frequencies, with a total emission of 
at least 2 x 10^^ erg (for Hq = 50 km s~^ Mpc~^ ) in 8.5 years. A model-dependent 
extrapolation of this value indicates that the total radiated energy may be as high as 
10^^ erg. The high value of the radiated energy supports a scenario in which most of 
the kinetic energy of the ejecta is thermalized and radiated in a short interaction with 
a dense circumstellar medium of nearly constant density. In this sense, 1988Z is not a 
supernova but a young and compact supernova remnant. 

Key words: supernovae: individual: SN 1988Z - supernova remnants 



1. INTRODUCTION 

The explosion of massive stars, thought to take place as the 
result of pair formation, collapse of a Fe core or neutron- 
ization of an ONeMg degenerate core, leads to supernovae 
(SN) that show H signatures in their spectra and thus are 
classified as type II. These objects display a wide range of 
observational properties which justifies the existence of the 
two photometric sub-classes, II-P and II-L (Barbon, Ciatti 
& Rosino 1979) , and the recent introduction of the peculiar 
spectroscopic sub-group Iln (Schlegel 1990). 

The spectra of SN Tin are characterized by the pres- 
ence of prominent narrow emission lines sitting on top of 
broad components of FWHM~ 15000 km/s at maximum 
light which look very similar to those of Seyfert 1 nuclei and 
QSOs (Filippenko 1989). They don't show the characteristic 
broad P-Cygni signatures of standard SN, although narrow 
P-Cygni profiles are detected in some cases at high spectral 
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resolution (SN 1997ab, Salamanca et al. 1998; SN 1995G, 
ESO data archive; SN 1997eg, Tenorio-Tagle G., priv. com- 
munication). SN Iln are normally associated with regions of 
recent star formation. Despite these general characteristics, 
SN Iln as a group exhibit some heterogeneity (see Filippenko 
1997, Turatto et al. 1999). 

SN 1988Z, one of the most extensively observed SN 
of this class, was discovered in 1988 December 12 in the 
galaxy M-l-03-28-022 (Zw 095-049) (Cappellaro & Turatto 
1988, PoUas 1988a). Earlier observations in late spring 1988 
show no evidence of the event (Turatto et al. 1993a). This 
is an exceptionally bright and peculiar SN in its spectro- 
photometric properties: 

• It is characterized by an extremely slow decay of lumi- 
nosity after maximum light, which makes it at day 600 ap- 
proximately 5 mag brighter in the V-band than standard 
SN II-P or SN II-L (Stathakis & Sadler 1991). 

• It has a strong Ha emission, with peak luminosities of 
about 4 X 10*^ erg s"^ (for Ho = 50 km s"^ Mpc"^ ) at day 
200 (Turatto et al. 1993a). This prodigious luminosity is 5 
orders of magnitude larger than that of SN 1987A. 

• Very high-ionization coronal lines (e.g. [Fe XJA6375 A, 
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[Fe XI]A7889-7892 A) are identified in the optical spectra, 
at least until day 492 (Turatto et al. 1993a). 

• At 2 to 20cm it is one of the most powerful radio-SN 
in the sky, with peak-luminosities up to 3000 times that of 
remnants like Cas A (Van Dyk et al. 1993). 

• Even 6 yr after maximum the SN shows an X-ray emis- 
sion of more than 10^^ erg s~^ (Fabian & Terlevich 1996). 

The summary of these points is that this is one of the most 
energetic radiative events known to be generated in a single 
stellar object. 

These radiative properties have been interpreted in the 
light of quick re-processing of the mechanical energy of 
a SN explosion by a dense circumstellar medium (CSM) 
(Chugai 1991, Terlevich et al. 1992, Terlevich 1994, Chugai 
& Danziger 1994, Plewa 1995). As explained in the ap- 
pendix, radiative cooling is expected to become important 
well before the thermalization of the ejecta is complete. As 
a result, the shocked material undergoes a rapid conden- 
sation behind both the leading and reverse shocks. These 
high-density thin shells, the freely expanding ejecta and the 
still unperturbed interstellar gas are all ionized by the ra- 
diation produced in the shocks, and are responsible for the 
complex emission line structure observed in these objects. 

In the case of SN 1988Z, a direct measurement of the 
CSM electron density was possible. The value determined 
from the [OIII]A4363 A to [OIII]A5007 A forbidden line ra- 
tio in the early stages of the evolution is between 4 x 10® 
and 1.6 x lO'^ cm"^ (Stathakis & Sadler 1991). This high 
density points towards SN 1988Z in particular, and SN Iln 
in general, being young compact SN remnants (cSNR). 

In this paper we present new optical and X-ray spectro- 
photometric observations of the late evolution (4.3 to 8.5 
yr) of this object, which combined with already published 
information on its early evolution, aims at deriving the to- 
tal energy radiated since discovery, and shedding new light 
into the mechanism that generates such a long-lasting bright 
event. 



2. DATA SET 

2.1. Optical observations 



The new optical observations that we present in this pa- 
per span a period of about 4 years, starting 4.3 years after 
discovery. 

B, V and i?-band photometry of SN 1988Z was ob- 
tained at the European Southern Observatory (ESO) in La 
Silla on six nights between 1993 March 27 and 1997 Febru- 
ary 10 using the 3.6m telescope in combination with the 
ESO Faint Object Spectrograph and Camera (EFOSCl) and 
the 3.5m New Technology Telescope (NTT) with the ESO 
Multi Mode Instrument (EMMI) or the Superb Seeing Im- 
ager (SUSI). 

When the nights were photometric the magnitudes were 
calibrated using a number of Landolt's (1992) standard 
stars. In non-photometric nights the SN magnitudes were 
derived using the local sequence defined in Turatto et al. 
(1993a). The SN magnitudes were measured by fitting a 



point spread function with the ROMAFOT package of MI- 
DAS. This technique allows a good subtraction of the back- 
ground and gives reliable results even for faint objects (see 
Turatto et al. 1993b). The journal of observations together 
with the resulting magnitudes and an estimate of the inter- 
nal errors are reported in table |^, where age refers to the 
time elapsed since the SN discovery, 1988 December 12 (JD 
24447508). 

The journal of the spectroscopic observations is given in 
table ^ where for each spectrum column 1 gives the date of 
acquisition, column 2 the age, column 3 the instrumentation 
used, column 4 the exposure time, column 5 the wavelength 
range and column 6 the resolution derived from the average 
FWHM of the night-sky lines. In two epochs we obtained 
spectra during two consecutive nights. In these cases we co- 
added the spectra to reach higher signal-to-noise ratios and 
included the cumulative exposure times in column 4. The 
spectra were wavelength calibrated with He-Ar arcs and flux 
calibrated using standard stars from the list of Hamuy et al. 
(1994). The absolute flux calibration of the spectra was ver- 
ified against the broad-band photometry and the agreement 
was generally fair. 

Figure 1 shows as an example of its late stages in evo- 
lution, the spectrum of SN 1988Z on 1996 February 14. The 
earlier evolution of the spectrum can be seen in Stathakis & 
Sadler (1991), Filippenko (1991) and Turatto et al. (1993a). 

2.2. X-ray data 

SN 1988Z has been observed three times with the 
Rontgensatellit (ROSAT) High Resolution Imager (HRI), in 
1995 (detection reported in Fabian & Terlevich 1996), 1996 
and 1997. Table 3 summarizes these observations. 

Only the first and last observations provide a good de- 
tection of the object; the middle one was too short. The 
count rate declined from 1.1 x 10~^ to 3.7 x 10~* s"^ be- 
tween 1995 May and 1997 May. Using Poisson statistics on 
the numbers of counts detected, we find that the probability 
that the average count rate was equal to or higher than that 
seen in 1995, and equal to or lower than that seen in 1997 is 
less than 0.2 per cent, adopting the weighted mean of both 
observations at both times. This indicates that the decrease 
is statistically significant. 

The 0.2-2 keV fiuxes, corrected for absorption in our 
Galaxy Nh = 2.5 x 10^° cm'^ (Burstein & Heiles 1982), are 
given in table 3, assuming a 1 keV bremsstrahlung spectrum. 
The 0.2-2 keV luminosity therefore dropped from (8 ±3.2) x 
10*" erg s"^ in 1995 to (3 ± 1.6) x 10*" erg s^^ in 1997. 

These are probably underestimates of the X-ray lumi- 
nosity. The hydrogen column density can be much higher 
than that assumed if the pre-ionized CSM is only partially 
ionized. Using the photoionization code CLOUDY 90 (Fer- 
land 1997) and a description of the evolution of the leading 
shock of SN 1988Z based on the model by Terlevich (1994) 
(see appendix) we estimate that if the mass of the CSM still 
to be swept is of the order of IM0, then the column den- 
sity would be of the order of 10^^ cm~^ at t — 12tsg, where 
tsg refers to the time of the radiative onset of the remnant 
(eq. A3), around 250 days from core collapse for this object. 
For a column density of lO^'^ cm~^ the previous estimates 
of the X-ray luminosity should be increased by a factor of 
eight if we assume a 1 keV temperature. The bolometric lu- 
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Table 1. Broad-band photometry 



date 


age 


B 


V 


R 


instrument 




days 


mag 


mag 


mag 






1993 Mar 27 


1565 


22.23 ±0.26 


22.08 ±0.25 


20.79 ±0.20 


NTT - 


h EMMI 


1994 Jan 14 


1856 






21.36 ±0.20 


3.6m ^ 


- EFOSCl 


1994 Dec 30 


2209 






21.71 ±0.20 


NTT - 


h EMMI 


1995 Jan 10 


2220 




22.73 ±0.20 


21.70 ±0.15 


NTT - 


h SUSI 


1996 Feb 14 


2620 






21.87 ±0.20 


3.6m ^ 


- EFOSCl 


1997 Feb 10 


2982 






22.04 ±0.20 


3.6m ^ 


- EFOSCl 



Table 2. Journal of spectroscopic observations 



date 


age 


instrument 


cxp. 


range 


rcsol. 




days 






min 


A 


A 


1993 Mar 27 


1565 


NTT - 


1- EMMI 


105 


3750-9700 


13 


1994 Jan 13-14 


1859 


3.6m H 


- EFOSCl 


115 


3700-6900 


16 


1994 Dec 30 


2209 


NTT - 


1- EMMI 


90 


3700-8900 


8 


1996 Feb 14 


2620 


3.6m H 


- EFOSCl 


120 


3750-6950 


16 


1997 Feb 9-10 


2982 


3.6m H 


- EFOSCl 


240 


3730-6900 


16 



minosities for the assumed spectrum are 1.6 times greater. 
We should emphasize that the absorption history is very de- 
pendent in the total CSM mass. For larger CSM masses the 
recombination occurs earlier in the evolution and thus the 
soft X-ray correction is larger. 

ROSAT observations of SN 1986J, another SN Iln, also 
show a high column density of lO^^cm"^ and a prominent 
X-ray luminosity 2 — 3 x 10*" erg s~^ about 9 years after the 
SN event (Bregman & Pildis 1992). 



2.3. Catalogue of radio, optical and X-ray 
photometry 

A compilation of photometric data gathered from the liter- 
ature, combined with the new spectro-photometric observa- 
tions introduced in sections 2.1 and 2.2 is presented in ta- 
ble 4. Column 1 gives the age since discovery. Columns 2 to 
10 give the flux and apparent magnitude of the object in the 
0.2 to 2 keV band from ROSAT; B, V, _R-band photometry 
and Ho flux from spectroscopy obtained with ground-based 
optical telescopes; and 2 to 20 cm fluxes from VLA. All 
fluxes are corrected for Galactic extinction using Seaton's 
law (1979). Column 11 gives the reference from which the 
data was extracted, where SS91 refers to Stathakis & Sadler 
(1991), T93 to Turatto et al. (1993a), FT96 to Fabian & Ter- 
levich (1996) and lAUC to the lAU circulars numbers 4691, 
4696, 4742 and 4761. (Cappellaro & Turatto 1988; PoUas 
1988a, 1988b, 1989; Gaskell & Koratkar 1989) 

Figure 2 shows a combined representation of the light 
curves of SN 1988Z at different frequencies. In the area oc- 
cupied by the optical and X-ray data, solid crosses corre- 
spond to Ha, solid squares to i?-band, solid triangles to V- 
band, open triangles to B-band and the solid hexagons to 
the ROSAT 0.2 to 2 keV band. 

The thin solid line that passes through the broad-band 
optical and X-ray points describes a t~^^^'^ law. This is a 
semi-analytic approximation of the bolometric emission pro- 
duced by a cSNR when it evolves in a dense CSM of homo- 



geneus density (eq. A7) having the canonical value of den- 
sity of the CSM 10^ cm~^, and being scaled to match the 
V-band luminosity of SN 1988Z. The overall decline is re- 
markably well followed by this law, even at time-scales for 
which the assumptions of the semi-analytic approximations 
are no longer valid {t ^ 250 days). The dashed line that 
passes through the Ho points represents the Ho luminosity 
predicted by the same models (Terlevich 1994), with no in- 
trinsic scaling. Although the models predict the peak of Ho 
luminosity correctly at day ~ 200, in the late evolution of 
the SN they fail to reproduce the more pronounced decay 
by about an order of magnitude. 

The lines that pass through the radio points are the 
models of an external thermal-absorbing gas combined with 
an internal thermal-absorbing/non-thermal-emitting gas fit- 
ted by Van Dyk et al. (1993) to their VLA observations. 

Since Ho, at the recession velocity of SN 1988Z Vr — 
6670 km s"^ (Barbon, Cappellaro & Turatto 1989), is fully 
included in the R pass-band we have estimated a Ho line- 
free _R-band magnitude subtracting an interpolated Ho flux 
from table 4 as 

i?co.^-2.51og(l0-/--^^^) , (1) 

where is the photometric zero of J?-band in the same units 
as f(Ho), and the 0.80 factor represents the transmission 
value of the redshifted Ho line relative to the peak trans- 
mission value of the _R-band filter used in the observations 
(see figure 3). The corrected _R-band magnitudes (-Rcont) are 
represented with empty squares in figure 2. 



3. SPECTRAL LINE ANALYSIS 

Figure 4 and 5 show the evolution of the Ho and H/3 profiles 
as a function of age, from days 155 to 2981. 

Following the analysis performed by Turatto et al. 
(1993a), we report in table ra the main emission line pa- 
rameters as measured with the ALICE package in MIDAS. 
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Table 3. ROSAT HRI data for SN 1988Z. The 5th column tabulates the probability that the counts seen 
arc a random fluctuation in the background. The flux in the 6th column has been corrected for Galactic 
absorption. 



date 


age 


count rate 


exp. 


random chance 


flux (0.2-2 keV) 




days 


10-4 s-i 


s 




X 10-^4 erg s-l cm--^ 


1995 May 16-25 


2355 


11 ±4 


12287 


2 X 10-5 


4± 1.6 


1996 Dec 14 


2924 


5±4 


6739 


0.1 


2± 1.6 


1997 May 13-24 


3085 


3.7 ± 2 


34300 


0.02 


1.5 ±0.8 



o ^ 

X 




4000 



5000 



6000 



7000 



A(A) 



Figure 1. Spectrum of SN 1998Z on 1996 February 14. 



This package allows multiple gaussian fitting of complex 
line profiles. The coding for broad, intermediate and narrow 
components follows that adopted in Turatto et al. (1993a). 
These were the components identified in the early evolution 
of SN 1988Z which we are following at later epochs. The 
analysis presented here is thus fully complementary with the 
results presented in that report. Figure 6 shows the decom- 
position performed in the Ho line of the spectrum taken on 
1996 February 14 as an example. In table pi a colon indicates 
that the value is uncertain. 

The decrease of the width of the broad component of 
Ho, shown in figure 7 with filled squares, is monotonical. 
The behaviour can be described by a power-law, which is 
similar to the time evolution of the velocity of the outer 
shock generated by a supernova remnant which evolves in a 
high-density medium, t'^^"^ (eq. A5), and which is plotted as 
a dashed- line. The law has been scaled to the velocity width 
of the line ~ 4000 km/s at day 1000. By comparison, the 
width of the intermediate component of Ho (plotted with 
empty squares) decreases in a much slower way. 

Figures 8 and 9 show the observed trend of the Balmer 
decrement and the [O III]A5007 A luminosity with time. 
Both show an initial increase followed by a sharp down turn 
and a slower decay. Also shown in Figure 8 are the predic- 
tions of the simple model of cSNR evolution. In qualitative 
terms the observations follow the trend prescribed by the 
simple model. The discrepancies, though, strongly suggests 
that there are fundamental aspects that are not considered 



by this model. These discrepancies together with the evolu- 
tion of the line profiles an multifrequency light curve, will 
be the basis for the construction of more detailed numeri- 
cal models which will include the role played by the density 
profile of the ejecta and the time dependent aspects of the 
CSM evolution. 



4. PHOTOMETRIC ANALYSIS 

4.1. Spectral Energy Distribution 

In figure 2 we observe that the light evolution in different 
bands peaks at different times: first in the optical, before dis- 
covery; then in Ho (i.e. in the ionizing continuum), between 
100 and 500 days; and lastly at radio-frequencies between 
400 and 2000 days. The spectral energy distribution (SED) 
thus changes with time. 

In order to include the ionizing radiation, which is ab- 
sorbed by the cool material, we estimate its value from the 
flux of a prominent recombination line. Ho is particularly 
useful, since several percent of the ionizing energy is re- 
radiated through this line. Assuming case B recombination 
at 10000 K, which is the average temperature of the gas that 
emits broad lines in the hydrodynamical models of cSNR 
(Terlevich et al. 1992), we can convert the Ho fiux to ion- 
izing flux. The ionizing flux /ion is related to the flux of 
ionizing photons ?(H'') by 
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Table 4. Photometric evolution from X-ray to radio 



Age /o.2 - 2 keV B V R /jja /2 cm /a.e cm fa cm /20 cm reference 

day erg cm^^ s^^ mag mag mag erg cm~^ mJy mJy mJy mjy 

16.54 T93 

2 16.59 16.23 15.93 T93 

4 16.59 15.98 T93+IAUC 

23 16.71 16.46 16.04 SS91 

24 16.69 T93 
34 1478x10-1'^ SS91 
50 16.84 T93 
53 17.04 16.79 T93 
57 16.68 16.35 lAUC 
59 16.82 16.37 lAUC 
62 1669xlO~i'5 SS91 
85 17.44 T93 
87 18.65 17.53 T93 
91 17.59 lAUC 

114 17.93 17.87 T93 

115 17.74 17.47 1449x10-1'' T93 
117 17.91 T93 
137 18.03 17.74 17.53 T93 
199 1630x10"" SS91 
282 19.02 17.80 T93 
392 0.67 vD93 
407 19.34 18.05 T93 
441 19.47 19.12 18.10 T93 
446 0.77 vD93 
451 19.19 T93 
464 1656x10-1'^ SS91 
474 18.09 T93 
492 1947x10-1"^ T93 
552 1.21 vD91 
557 19.54 18.38 SS91 
596 1.26 2.10 vD93 
653 1.38 1.15 vD93 
684 1019x10-18 T93 
717 936x10-1'^ T93 
731 20.78 20.16 T93 
736 20.25 T93 
738 730x10-1'^ T93 
743 20.79 20.09 18.88 T93 
751 1.78 2.09 1.43 vD93 

758 793x10-16 T93 

759 18.82 T93 
798 20.89 20.24 18.94 814x10-1'^ T93 
800 19.11 T93 
845 21.29 20.61 19.08 T93 
930 0.47 1.90 1.68 1.31 vD93 

1028 1.85 1.48 0.59 vD93 

1086 20.82 19.53 T93 

1095 20.59 19.48 T93 

1148 21.47 20.89 19.76 T93 

1149 338x10-1*' T93 
1160 1.05 1.67 1.50 1.17 vD93 
1420 0.91 1.57 1.22 0.91 vD93 
1487 1.45 1.72 1.04 0.78 vD93 
1565 22.02 21.92 20.67 124x10-1*^ this paper 
1856 21.24 84.7x10-1" this paper 
2209 21.59 56.0x10-1'' this paper 
2220 22.57 21.58 this paper 
2355 4x10-1^^^ FT96 
2620 21.87 37.6x10-1" this paper 
2924 2x10-1'* this paper 
2982 22.04 22.2x10-1" this paper 
3085 LSxlO-l* this paper 
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-18 I ^-^ ^ ' ^ ' ^ 1 

1000 2000 3000 

age (days) 

Figure 2. Evolution of the ufi^ light curve in radio to X-ray bands. In the upper part of the panel open triangles correspond to _B-band, 
solid triangles to V-band, solid squares to i?-band, empty squares to i?cont band, solid hexagons to the ROSAT 0.2 to 2 keV band and 
crosses to Ha. The thin solid and dashed lines are models for the bolometric light (scaled to the values of the optical data) and Ha 
evolution of a cSNR. In the lower part of the panel, the VLA radio data at 2, 5.5, 6 and 20cm is represented with empty squares, empty 
triangles, solid squares and empty triangles. The line fittings correspond to the models of Van Dyk et al. (1993). The thick solid line on 
the top of the panel shows the strong lower limit estimate while the thick dashed line represents our best estimate of the bolometric light 
curve of SN 1988Z (see section 4.2). 



Table 5. Spectral line decomposition of the late spectra of SN 1988Z 



age 






Ho 




H/3 


Hel-NalD 


days 




broad 


interm. 


interm. 


narrow 




1565 


FWHM (A) 


48.6 


25.1 






75: 




flux (xlO-l^ ergs-lcm-2) 


26 


98 






6: 


1859 


FWHM (A) 


40.7 


28.3 


60.9 


16.3 


25.3 




flux (xlO~^'' ergs~^cm~-^) 


8 


67 


2.7 


3.0 


2.2 


2209 


FWHM (A) 


49.2 


27.3 




18.6 


47: 




flux (xlO-l^ ergs-lcm-2) 


5 


47 




2 


1.4: 


2620 


FWHM (A) 


45.1 


28.6 


61.1 


16 


21: 




flux (xlO-i^ ergs-icm-2) 


9.6 


24 


1.7 


1.9 


0.4: 


2981 


FWHM (A) 


53.9 


25.3 


45.7 


16.7 


51.2 




flux (xlO-i^ ergs-icm-2) 


6.9 


12.9 


1.2 


1.3 


1.0 
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500 600 700 800 900 1000 

A(nm) 

Figure 3. Transmission values of tlic three _R-band filters used in the observations (ESO#554, ESO#608, ESO#642). Also plotted is 
the Ha spectral region of SN 1988Z on 1989 April 5 
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Figure 4. Evolution of Ha profile 
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Figure 5. Evolution of H/3 profile 



/ion = q{'H°)hi/u 



(2) 



where TTkm is the mean frequency of the ionizing photons, 
approximately 27 oV. For case B recombination, the Ha 
flux is related to the flux of ionizing photons by /jj^ = 
1.367088 X 10"^^g(H°) in c.g.s. units (Kennicutt 1988, Os- 
terbrock 1989). Thus the conversion between the fluxes is 



/ion « 32/jj^ 



(3) 



Figure 10 shows the SED at six stages in the evolution 
of SN 1988Z. Solid squares correspond to the flux values 
which were interpolated from those in table 4, and in the 
case of the ionizing continuum, estimated using equation 3. 
Open squares indicate that we have performed an extrap- 
olation using the fits of Van Dyk ct al. (1993) for radio 
wavelengths, or a f"^^'''^ law drawn from the nearest value 
of optical points and from the Ha or X-rays for the hard 
radiation. Therefore, open symbols should be taken with 
caution. The only exception being values extrapolated from 
data which arc ±50 days away, in which case these arc rep- 
resented as solid squares too. Crosses correspond to the Ha 
line-free Jicont flux. 

The SED of SN 1988Z differs from the SED of classical 



SNe in that its high energy distribution is very prominent. 
There is also an indication that the shape in the early evo- 
lution differs from the general assumption that the emer- 
gent spectrum is created by an expanding hot atmosphere 
with mainly a Plamckian output (e.g. Catchpole et al. 1989, 
Eastman & Kirshner 1989). The dashed line in Figure 10 
shows what the SED would look lino if the optical colours 
of SN 1988Z were to be interpreted as a blackbody distribu- 
tion, as has been assumed in the case in other type Iln SNe 
(e.g. Wegner & Swanson 1996). Hard-radiation and radio 
excesses are expected when the SN shock interacts with a 
CSM. At this stage we are viewing the kinetical reprocessing 
of a SN remnant 



4.2. Radiated energy estimates 

In order to calculate the energy radiated by SN 1988Z in the 
first eight years of its evolution, first we integrate the light 
curves of Figure 2, interpolating linearly in the time-axis be- 
tween the observed points, i.e. with no further extrapolation 
on the behavior of the light curves after the last observa^ 
tional point. This will give us an estimate of the observed 
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Figure 6. Multicomponent fit of the Ha line of the spectrum taken on 1996 February 14. The parameters of the components are given 
in tabic 5. The dot-dashed lines represent single gaussian components, while the solid line and dotted line represent the spectrum and 
the fit. 
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Figure 7. Evolution of the Ha line profile. Solid squares represent the broad component and empty squares represent the intermediate 
component. The dashed line represents the law for the velocity shock evolution of a cSNR. 



radiated energy during those eight years, which should be 
regarded as a lower limit to the total energy radiated during 
that time. Table 6 gathers the results, where time refers to 
the maximum span of the light curve in that band. Prom 
this table we derive that the observed emitted energy of 
SN 1988Z is dominated by the optical to X-ray radiation, 
having emitted at least 1.0 x 10^^ erg in eight years. 

This is a very rough estimate, which leaves gaps in the 
spectral energy distribution and in the time evolution at the 



different bands. As an example, if we interpolate linearly the 
monocromatic flux at frequencies between B and V and V 
and J?-bands along the evolution of the light curve, i.e. 

/■8.2yr i-i/^ 

Sopt = 4.n(f / dt dvU , (4) 
where d is the distance to the cSNR, then we obtain that 
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Figure 9. Evolution of [O III]A5007 A line 



the energy radiated from B to R band is 1.6 x 10 erg, or 
1.5 X 10®" erg if using -Rcont instead of R. 

In order to obtain a more complete estimate of the emis- 
sion at optical-UV wavelengths, we have performed an ex- 
trapolation of the B, V and i?cont fluxes to cover the 912A 
to 1/^m range, assuming that the SED is a power-law, i.e. 

(X v", where a is calculated from the B — V colours for 
the UV to F-band interval (typically —3.0 < a < —0.5) and 
from the V — i?cont colours for the V-band to l/itm interval 
(typically —2 < a < 5). If the color doesn't change dramati- 
cally between these regions, the optical energy inferred is of 
the order of 3.2 x 10^° erg for the eight year evolution. 

The ionization energy, in turn, can be estimated from 
the integration of the Ha light curve. The use of eq. 3 gives 
a value of 8.7 x 10™ erg. 

To estimate the radiation at 0.2-2 keV, we have extrap- 
olated the three X-ray points to earlier epochs than 6.5 yr, 
when the first ROSAT observation was obtained. If we adopt 



the t~^^^'' law before day 2355, so that at day 2355 it repro- 
duces the first ROSAT observation, wc have that the energy 
radiated in 3000 days in the 0.2-2 keV band could be as high 
as 6.0 X 10®° erg. If a complete bremsstrahlung spectrum of 1 
keV is considered, the energy emitted in X-rays alone could 
be about 9.6 x 10®° erg, and if we consider internal absorp- 
tion of the GSM, this could go up to 7.6 x 10®^ . 

These estimated values of the emitted energy in the 
8.5 yr of evolution of SN 1988Z at different wavelengths are 
summarized in Table 7. Beware that some entries in this 
table are subsets of other entries (e.g. 4000-7000A and 912- 
l/^m). and that some entries are in part the result of re- 
processing the energies of other entries (e.g. 912-1/im and 
bremsstrahlung). Considering that half of the leading shock 
radiated energy is emitted outwards and the other half in- 
wards and is therefore reprocessed by the high density thin 
shell, the total energy radiated in the event should include 
either twice the observed X-ray luminosity or the sum of 
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Table 6. Integrated values of observed emitted energy in different 
wavelengths 



band 


time 

yr 


energy 
erg 


0.2-2 keV 


2.0 


3.6 


X 


1048 


ionization 


8.2 


8.7 


X 


1050 


B 


4.3 


4.0 


X 


1049 


V 


6.1 


2.6 


X 


1049 


R 


8.2 


5.8 


X 


1049 


-Rcont 


8.2 


3.7 


X 


1049 


2 cm 


2.3 


2.4 


X 


1046 


3.6 cm 


2.4 


2.2 


X 


1046 


6 cm 


3.0 


1.4 


X 


1046 


20 cm 


2.3 


1.4 


X 


1045 



Table 7. Integrated values of estimated emitted energy in differ- 
ent wavelengths during the 8.5 yr of evolution 



band 


cnorg 

t'lg 


y 




bremsstrahlung 


9.6 X 


1050 


to 7.6 X 10'''^ 


0.2-2 kcV 


6.0 X 


105O 


to 4.8 X 10^1 


ionization 


8.7 X 


1050 




4000-7000A 


1.5 X 


1050 




912A- 1/im 


3.2 X 


1050 





the observed X-ray luminosity and the inferred ionizing and 
optical luminosity, or twice the sum of the observed ionizing 
and optical luminosity. 

A strong lower limit of the total radiated energy is then 
2 X 8.7 X 10^° -I- 3.2 X 10^° « 2 x lO'^^ erg. The bolometric 
light curve derived for this estimate is represented with a 
thick solid line in figure 2. 

A best estimate of the total radiated energy requires 
a model-dependent extrapolation either from the X-ray or 
from the Ha radiation. The X-ray derived values axe rather 
uncertain, since they require the assumption of time, spec- 
tral and absorption corrections. Allowing for these uncer- 
tainties, the energy could be as high as 2 x 7.6 x lO^'^ « 
10^^ erg. The bolometric light curve derived for this model- 
dependent estimate is represented with a thick dashed line 
in figure 2. Note that the excess energy of this second esti- 
mates is produced mainly at early epochs. A high value of 
the X-ray radiated energy is supported by the ionization en- 
ergy estimates. The models indicate that during the earlier 
stages of evolution the emitted radiation is very hard and 
most of it leaves the remnant without producing ionization. 
A rough estimate based on the simple model indicates that 
the total energy radiated is about six times larger than the 
ionizing radiation, that is about 4 x 10®^ . Recalling that this 
includes only half of the total energy emitted, we conclude 
again that the total energy radiated may be close to lO^'^. 

We emphasize that these estimates totally disregard any 
intrinsic reddening produced in the host galaxy or in the 
CSM of SN 1988Z. 



5. DISCUSSION 

We have shown that the integrated light emitted in radio, 
optical and X-rays by SN 1988Z in the 8.5 yr of evolution 
after discovery is of the order of several times 10^^ erg and 
may be close to 10^^ erg, i.e. larger than the canonical value 
of the kinetic energy released in a SN explosion, while typi- 
cal SNe radiate one to two orders of magnitude less energy 
(e.g. Panagia et al. 1980, Catchpole et al. 1987, 1989). The 
value of the radiated energy in SN 1988Z is slightly smaller 
than that estimated for the most energetic kinetical release 
in a SN explosion observed so far, (2 — 5) x lO'"''^ erg, in 
SN 1998bw (Galama et al. 1998, Iwamoto et al. 1998). In- 
deed on the basis of our result we speculate that SN 1988Z 
was a hypernova and may have been caused by the formar 
tion of a black hole in the core of a massive star. Most of 
the energy in SN 1988Z is radiated in the optical to X-ray 
interval, with spectra, SED and light curves departing from 
those of classical SN II. The high densities derived from the 
O forbidden lines and the flat light-curves strongly support 
the idea of a high-density circumstellar shell reprocessing 
the kinetic energy through fast radiative shocks. 

The estimated CSM mass swept by the shock in this 
time is of the order of 15Mq (eq. A8), although this quantity 
is dependent on the assumed model. Chugai & Danziger 
(1994), for example, give an estimate of only about IMq 
from their own models. 

While a detailed comparison of light-curve, SED, and 
line profile evolution is still missing, we have shown that 
the simplified model of ejecta interacting with a dense and 
homogeneous CSM outlined in the appendix reproduces the 
overall decay of the optical light, the level of X-ray and Ha 
emission and its line-width evolution. The evolution of the 
Ha flux is however more rapid than that predicted by the 
canonical models, and the absolute values of line- widths and 
Balmer decrements differ sometimes by more than a factor 
of 2. More fine-tuned models, possibly introducing a small 
density gradient in the circumstellar shell, arc thus required. 

Observatioiially much work is also needed, on one hand 
to identify cSNR like SN 1988Z early enough in their evo- 
lution, and on the other to obtain good quality data over 
the widest possible spectral range and time coverage. High- 
quality and high-resolution spectral information is particu- 
larly needed in the X-ray region early enough in the evo- 
lution of the SN, where and when most of the energy is 
radiated. A close optical monitoring of these objects would 
allow us to assess the predicted existence of modulations 
in the light curve decay and line component anomalities in 
the spectra due to catastrophic radiative cooling, and the 
precise dating of the disappearance of the broad emission 
lines. 
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APPENDIX A: THE "SIMPLE MODEL" FOR 
COMPACT SUPERNOVA REMNANTS 

The interaction of the SN ejecta with a dense CSM causes 
a shocked region of hot gas enclosed by two shock waves: on 
the outside the leading shock, and on the inside the inward 
facing "reverse" shock. The leading shock {vg ~ iC km s~^) 
encounters dense circumstellar material and raises its tem- 
perature to ~ 10® K. The reverse shock, which is initially 
substantially slower {v ~ 10^ km s~^), begins to thermal- 
ize the supernova ejecta to temperatures of about lO'^ K. 
Early analytical and numerical computations of the evo- 
lution of SNRs in a dense medium (Chevalier 1974; Shull 
1980; Wheeler et al. 1980) showed a speeded-up evolution 
compared with the "standard" solution in a medium of no 
= 1 cm~^. All evolutionary phases (free expansion, thor- 
malization of the ejecta, the quasi-adiabatic Sedov phase, 
the radiative and the pressure modified snow-plough phases) 
which have been thoroughly studied for the standard case, 
are substantially speeded up as a consequence. 

Supernova remnants evolving in a dense CSM reach 
their maximum luminosity (L > 10* L©) at small radii 
(-R < 0.1 pc ) soon after the SN explosion {t < 20 yr) while 
still expanding at velocities of more than 1000 km s"^ (Shull 
1980; Wheeler et al. 1980; Draine and Woods 1991; Ter- 
levich et al. 1992; Terlevich 1994). In these cSNR, radiative 
cooling becomes important well before the thermalization 
of the ejecta is complete. As a result, the Sedov phase is 
avoided and the remnant goes directly from the free expan- 
sion phase to the radiative phase. The shocked matter un- 
dergoes a rapid condensation behind both the leading and 
the reverse shocks. As a consequence two high-density, fast- 
moving thin shells are formed. These dense shells, the freely 
expanding ejecta and a section of the still dynamically un- 
perturbed interstellar gas, are all ionized by the radiation 
from the shocks. 

For SNe evolving in CSM densities of the order of 
no ^ 10^ cm~^, the cooling-time (tc) and cooling-length 
(tc) scales for the post-shock temperatures are very small. 



3kTs . ^ vj 

tc — T. r — 0.2 f — yr, 

8noA nyAas 

1 7/^ 

rc ~ -tcVs ~ 1.8 X 10" — I— cm, 
4 n7A23 



(Al) 



(A2) 



in 10^ cm~"^units, A23 is the cooling function (A) in lO^"^"* erg 
cm* s~^ units, Ts is the shock temperature, and k is Boltz- 
maim's constant. Thus, radiative losses become important 
at very early times when the shock velocities and tempera- 
tures arc Vs > 10'^ km s^^ and Ts > 10^ K, well before the 
ejecta is even thermalized. This means that a large flux of 
ionizing photons will emerge from the shocked gas at X-ray 
energies. 

For a supernova remnant which injects 10®^ ergs into a 
medium of constant density 117, the onset of the radiative 
phase behind the leading shock which causes the formation 
of a dense outer shell, assuming free-free cooling only, (Shull 
1980, Wheeler et al. 1980, Draine and Woods 1991) begins 
at a time tsg, given by 



tsg = 230 -Egj^ 



1/8 -3/4 



days 



(A3) 



where S51 is the energy deposited by the SN in units of 
10®^ ergs. At this stage, the shock is at a radius 



with velocity 



4600 E, 



1/8 1/4 



51 



-5/7 



km s 



temperature 

Ts = 3.0 X 10« E'l^ (—) K 

V tsg J 

and luminosity 



= 2 X 10« Eli^ r^^ 



-11/7 
— ) ergs 



The mass swept by the shock is 



M. 



,.e« = i.i£;f/%-v^(^) Me. 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



where is the velocity in 10® km s ^ units, n^ is the density 



These approximate formulae assume that the ejecta has 
already been fully thermalized. However, for the case of in- 
terest in this work (i.e. for no > 10® cm~*) strong radiative 
losses occur before thermalization is complete. Because the 
cooling processes radiate the thermal energy at the same 
rate as thermalization proceeds, the Sedov phase is totally 
inhibited and thus there is no self-consistent analytic treat- 
ment for the evolution of such remnants. It is therefore nec- 
essary to follow the detailed time evolution of the gas flow. 
In particular, special care should be given to the post-shock 
structure which is sensitive to the details of the ambient 
density distribution and to the temperature dependence of 
the radiative cooling function. 

The cooling time scale, tc, for am optically thin plasma, 
is proportional to the inverse of the gas density and the evo- 
lution proceeds faster at higher ambient densities. On the 
other hand, the temperature dependence of the coofing func- 
tion is different for different temperature ranges. Adiabatic 
shocks are stable but cooling instabilities can develop over 
a wide range of radiative shock conditions (Avedisova 1974, 
Falle 1975, 1981, McCray, Stein and Kafatos 1975, Chevaher 
and Imamura 1982, Imamura 1985, Bcrtschinger 1986). 

The details of the transition from a nearly adiabatic to 
a strongly radiative shock depend on the ability of the gas 
to readjust to the cooling rate. Pressure gradients tend to 
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be smoothed out in a sound-crossing time, td, and the ratio 
tc/td provides an estimate of the conditions prevailing in the 
cooling gas. For tc/td > 1, at moderate cooling rates, the gas 
elements are continuously compressed as their temperature 
falls and the cooling process operates quasi-isobarically at 
the pressure attained by the gas immediately behind the 
shock. For tc/ta < 1, however, the cooling rate dominates 
over any pressure readjustment and the process becomes 
quasi-isochoric at the post-shocked density of the cooling 
gas elements. 

A large pressure imbalance then develops in the flow, 
and new additional shocks are generated which end up com- 
pressing the cooled gas. This process, termed "catastrophic 
cooling" (Fallo 1975, 1981), appears during thin shell for- 
mation and the instabilities continue to operate during the 
rest of the radiative shock evolution (Chevalier and Imamura 
1982; Bertschinger 1986; Cioffi, McKee & Bertschinger 1988; 
Tcnorio-Tagle et al 1990). The catastrophic cooling acquires 
a central role in the case of supernovae evolving in high 
density media due to the strength of the radiation produced 
upon cooling, and the rapid variations inherent in the shock 
propagation. These features imply that a large flux of ion- 
izing photons will emerge from the shocked gas. The wide 
range of gas temperatures in the cooling region results in a 
power-law-like spectrum at UV and X-ray frequencies (Ter- 
levich et al. 1992). 
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